Non-classical photon streams using rephased amplified spontaneous emission 
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We present a fully quantum mechanical treatment of optically rephased photon echoes. These 
echoes exhibit noise due to amplified spontaneous emission, however this noise can be seen as a 
consequence of the entanglement between the atoms and the output light. With a rephasing pulse 
one can get an "echo" of the amplified spontaneous emission, leading to light with nonclassical 
correlations at points separated in time, which is of interest in the context of building wide bandwith 
quantum repeaters. We also suggest a wideband version of DLCZ protocol based on the same ideas. 

PACS numbers: 03.67.-a, 32.80.Qk, 42.50 p 



I. INTRODUCTION 

In order to extend the range of quantum key distribu- 
tion, quantum networks and tests of Bell inequalities, a 
method for efficiently generating entanglement over large 
distances is required. To achieve this goal a quantum 
repeater is necessary Such repeaters are generally 
based on methods for entangling one light field entangled 
with another at a later point in time. This has led to in- 
creasing interest in quantum memories for light, which in 
conjunction with pair sources would achieve this. Many 
impressive experiments have been performed in the area 
of quantum memories and repeaters. The quantum state 
of a light field has been stored in a vapour cell with high 
fidelity, and then measured at a later time • Single pho- 
tons and squeezed states have been stored and recalled 
S 0i Q , and non-classical interference of the light from 
distant ensembles has been observed @. Entanglement 
of and teleportationQ between two distant trapped 
ions has been achieved using an optical channel, using 
DLCZ type measurement induced entaglement. 

Photon echoes have a lon g hi story of use in classical sig- 
nal processing with light [qLIIOII . There are now a number 
of proposals and experiments [ll|, [H, [H, [H, [l^, [lB| re- 
lated to the development of photon echo based quantum 
memories. A distinct advantage of echo based techniques 
is that they arc multimode 17]. 

Current photon echo quantum memory techniques, all 
involve some modification of the inhomogeneous broad- 
ening profile. This imposes limits on the range of suit- 
able materials. It would be much more convenient to use 
something akin to the standard two pulse echo as a quan- 
tum memory which doesn't require such modification. 

The paper is arranged in three sections. In Sec. Hi] 
we present the quantum mechanical Maxwell-Bloch equa- 
tions for atomic and photonic fields. Then in Sec. IIIIl 



as an example of this formalism, we present an analysis 
of the standard two pulse photon echo and its applica- 
bility as a quantum memory. The two pulse echo as a 
quantum memory has already been investigated by oth- 
ers IS]. We revisit the problem here and show that this 
protocol fails as a quantum memory due to the strong 
rephasing pulse, which inverts the medium and causes 
additional noise on the output photonic fields. Finally 
in Sec. IIVI after exploring the origin of this noise, we 
propose that this noise can be rephased to lead to time 
separated, temporally multimode, wide bandwith photon 
streams with non-classical correlations. So while a stan- 
dard two pulse echo fails as a quantum memory, rephased 
amplified spontaneous emission (RASE) can be utilized 
in the DLCZ protocol JJ,]. With this modified DLCZ 
protocol, the inhomogeneous broadening no longer limits 
the time between the write and read pulses but instead 
increases the bandwidth of the process. This is of signif- 
icance to current experiments, were the inhomogeneous 
broadening is an issue [H [HI, [H, Hi] . 



II. QUANTIZED MAXWELL-BLOCH 
EQUATIONS 

We shall model an inhomogeneously broadened col- 
lection of two level atoms interacting with a 1-D field 
propagating in one direction, with the following quan- 
tum Maxwell-Bloch equations: 

d 

— a-iz,A,t) =zAo-_(z, A,t) - i a{z,t) a^iz, A,t) (1) 
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where (T+,_,2 represent the quantum atomic spin oper- 
ators, d is the quantum optical field operator, a is the 
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FIG. 1: (Color online) A two pulse photon echo sequence 
showing the approximations made in the treatment, a weak 
first pulse is applied to the system and is recalled using and 
ideal 7r-pulse. 



optical depth parameter, which depends on the coupling 
between the atoms and the field and on the atom den- 
sity. The parameter A is the detuning from some cho- 
sen resonant frequency and z is the distance along the 
propagation direction. The operators have the following 
commutation relations: 



[(7,(z,A,t),aj(z',A',i)] 
x6{z- z')S{A~ A') 



2tt . , , ^ 
— <^ijkcrk[z, A, t) 



(4) 



(5) 



As can be seen from Eq. [31 we take the density of atoms 
as a function of frequency to be a constant. In the 
case of rare earth ion dopants, the inhomogeneous broad- 
ening can be many times larger than the homogeneous 
linewidths, and as a result in most experiments without 
holcburnt features this is a good approximation. 

The above Maxwell-Bloch equations can be derived by 
dividing the atomic ensemble into thin slices and then 
modelling each slice as a small collection of atoms inside 
a Fabry-Perot cavity, using standard input-output theory 
[2^. Taking the limit as reflectivity of the mirrors go to 
zero one arrives at Eq. (IT]-[3]), where a{z,t) is the input 
field at the left hand side of the cavity and a{z + dz,t) 
is the output field at the right hand side of the cavity. 



d_ 

di 



Dg{z,A,t) = iADg{z,A,t)+ia{z,t) (6) 

d ia f°° - 

g-^a{z,t) = ^ Dgiz,A,t)dA. (7) 



Eq. [6]is just a first order linear equation with solution, 
Dg{z, t, A) = ~i f dt' a(z, t')e'^^*"*'^ 

J —00 

+ e'^'Dgoiz,A) , (8) 



where DgQ(z,A) is an initial condition. Taking the 
Fourier transform of Eqs. [7] and \S\ and substituting, one 
arrives at the following expression: 



—a{z,uj}^— dAa[z,uj) 
dz 2tt 



1 



+ nS{uj - A) 



V2^ J- 



i{uj - A) 
dAS{uj~A)DgQ{z,A) , 



—a{z,uj) + -j=Dgo{z,uj), 



(9) 



where S(uj) is the Dirac delta function. Solving Eq. [5] 
and fourier transforming back to the time domain we get 



d{z,t) = a(0,t)e~"^/2 



ta 

\/2tt Jo 



dz' e"(^'-^)/2 f)^^^,'^t) , 



(10) 



where a(0, t) denotes the input photonic field. Eqs.[8]and 
[TU]form the ground state solutions for all input times. 

After the 7r-pulse the atoms are all very close to the 
excited state {az ~ +1) in which case we can approxi- 
mate (7+ by a harmonic oscillator field D,,- This gives us 
the excited state Maxwell-Bloch equations. 



III. THE TWO PULSE PHOTON ECHO 

The first application of our quantum Maxwell-Bloch 
equations will be in analysing a memory based on a 
two pulse photon echo. The Maxwell-Bloch equations 
are non-linear and in general difficult to solve analyti- 
cally, however following work done with the semiclassical 
Maxwell-Bloch equations i26j one can make reasonable 
approximations that simplify the situation greatly. First 
we shall assume the input pulse is weak and is much 
smaller than a tt pulse. In this case all the atoms will 
stay near their ground state {az ~ —1) and we can ap- 
proximate the atomic lowering operator cr_ as a harmonic 
oscillator field Dg. The result are linear equations which 
we shall refer to as the ground state Maxwell-Bloch equa- 
tions. 



-Dl{z,A,t) = tADl{z,A,t)-id{z,t) (11) 
g^^M = — J^^Dl{z,A,t)dA. (12) 

We treat the tt pulse as being a perfect tt leading to the 
transformation <— Dg. We will discuss the treatment 
of the perfect tt pulse later in the text. 

Bringing Eqs. II II and ll2l through the same mathemati- 
cal process as Eqs.|6]and[71 we arrive at the excited state 
solutions: 

Dliz, t,A)^i f dt' a(z, t')e'^^'~''^ + e'^'bl^{z, A) , 

(13) 
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o(z, t) = a(0, t)e 
ia 



az/2 



'211 



dz'e"(^-^')/2 7)^o(z',t) , (14) 



where a(0, t) and d\q{z^ A) are initial conditions for the 
photonic and atomic excited fields respectively. 

Matching the ground (Eqs. [5] and [TU)) and excited 
(Eqs. [13] and [14]) state solutions at the point the tt- 
pulse is applied we get a complete solution. The ef- 
ficiency is sinh^(^) and in the limit of large optical 
depths high efhciencies are possible. Physically, this is 
because the photon echo is produced in the first piece 
of the sample and then gets amplified as it propagates 
through the inverted medium. The noise on the out- 
put can be quantified by considering the case of no in- 
put pulse, then the output will be amplified spontaneous 
emission (ASE), simply the vacuum noise amplified by 
the gain of exp(az) of the inverted ensemble. In the case 
of no input pulse, we get an incoherent output field with 
(at(t)a(t')) = <J(t - t')(exp(aO - 1). 

It is interesting to consider the source of this noise. In 
the model we have no dissipation and so the total system 
evolves through pure states. 

Eq. [5] is analogous the output of a beam splitter. The 
input fields being light and atoms, with the output fields 
consisting of combinations of photonic and atomic exci- 
tations. One can see that the addition of atomic excita- 
tions in the solution is necessary for the conservation of 
the commutation relations due to the input photonic field 
decaying away at large al. The excited state solution is 
analogous to a non-degenerate parametric amplifier [27| , 
here the input field is amplified, and the comutation re- 
lations are preserved by the addition of atomic creation 
operators. The state of one output mode is only mixed 
if the other is traced over, if the system is viewed as a 
whole one has an entangled state. In the next section we 
show that by applying a rephasing pulse to the ensemble 
we can turn the excitation of the atoms back into light, 
leading to streams of photons with highly non-classical 
correlations between two points separated in time. 
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inverts ions 
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Time 

"eclio" of amplified 
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FIG. 2: (Color online) Two tt pulse photon echo sequence pro- 
posed for generating rephased amplified spontaneous emission 
(RASE). 




FIG. 3: (Color online) Plot of R{al), show the violations of 
the Cauchy-Schwartz inequality for small optical depths. 



gion 2 we get from the final condition for region 1: 



Dgo(z,A) =ie""/2 / dt' a^O,t') 



iAt' 



— oo 




dt' e*^*' / dz' e"'^'-''^/^ Deo{z',t') 
27r 7_oo Jo 

(15) 



These boundary and initial conditions substituted in 
Eqs. ([S HU]) give a complete analytic solution of the lin- 
earised Maxwell-Bloch equations. 

To show that the photon streams described by these 
solutions have non-classical correlations we consider, 



IV. REPHASED AMPLIFIED SPONTANEOUS 
EMISSON 



Now we consider the two tt pulse sequence shown in 
Fig. [2 For region 1 the atoms will be inverted due to 
the first tt pulse and hence Eqs. pTl [T^ will apply. For 
region 2 the atoms will be near the ground state due to 
the refocusing tt pulse, hence Eqs. ^ [7]) describe the 
dynamics. We take the second 7r-pulsc to occur at i = 0. 

The solution for the light in region 1 is given by Eqs. [5] 
and [To] and the solution in region 2 is given by Eqs. [13] 
and 1141 For boundary conditions we take the incident 
field, a(0, t), to be in it's vacuum state as we do for the 
initial condition Deo(z, A). The initial condition for re- 



R: 



p{ti,ti)p{t2,t2)' 



(16) 



where p{ti,tj) = {a^l,ti) a{l,ti) cv {l,tj) a{l,tj)) . For 
classical fields the Cauchy-Schwartz inequality states 
that R < 1 [13 • Considering times equally separated 
about the second 7r-pulse, from the expression for the 
output fields derived above we get 



R{al) = 



al 



sh{al) 



4 sinh(f ) [e"' - 1] 



(17) 



Figure [3] shows that for small optical depths {al < 1) 
the output at times equally separated from the refocus- 
ing TT-pulse has nonclassical correlations. It should be 
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pointed out that the entanglement between times is not 
perfect because the echo efficiency is not 100%. The de- 
tection of a photon in Region 2 at a particular time means 
that there must have been one in Region 1 at the match- 
ing time, however the converse is not true. 

An advantage that RASE has is it's potential im- 
plementation in a larger range of systems. This is 
in contrast with the implementation of current pho- 
ton echo quantum memories. Current quantum mem- 
ory echo techniques use very fine spectral features pre- 
pared in the inhomogeneous line, rather than the natu- 
ral inhomogeneous profile and optical rephasing pulses. 
Indeed the fact that AFC[l5|] and CRIBJT,, 12.] type 
echoes ideally want homogeneously broadened ensembles, 
has led to their investigation in non-solid-state systems 
[1^. When selecting a rare earth ion system, one finds 
that the systems long lived spectral holes, such as eu- 
ropium or praseodymium, have inconvenient wavelengths 
(~ 580 nm and ~ 606 nm), in systems which are much 
more compatible with optical fibers and diode lasers the 
holes are much more transient, making high fidelity op- 
eration difficult at best. 



IMPERFECT n PULSES 



—Dg{z,p,A,t) = iADgiz,p,A,t) 



47r2 



—a{z,kt,t) 
oz 



//•OO 
d?p j dA Dg{z,p,A,t)e-^'^^r (19) 

Fourier transforming the atomic operators along the 
transverse dimensions by defining 



Dg{z,-kt,A,t)^ J d^pDg{z,p,A,t)exp{~i-kt.p) (20) 

leads to Maxwell Bloch equations which are diagonal in 
the transverse wave vector 



—Dg{z,kt,A,t) = iADg{z,kt,A,t) +ia{z,kt,1(]21) 

— a{z,kt,t) = — / dADg{z,kt^,t). (22) 
oz 2tt ,1^ 



So far we have treated the tt pulses as ideal and the 
effect of non-ideal tt pulses needs to be considered. It 
is feasible to make a pulse such that afterward one can 
make the approximation tr^ « 1 especially as we are in- 
terested in optically thin samples. The ability to do this 
in optically thick samples is also helped by the area the- 
orem, which states that a tt pulse remains a tt pulse as it 
propagates through a medium [3 ■ In the situation 
where tr^ « 1 after the pulse, we can model our non-ideal 
TT pulse as the combination of an ideal tt pulse and some 
excitation of the De field. This excitation of the field 
will be temporally brief, and if the inhomogeneous broad- 
ening is flat the ensemble of atoms will quickly dephase 
leading to no net polarization in the ensemble shortly 
after the n pulse. This means that the excitation pro- 
duced by the imperfect tt pulse will no longer interact 
with the optical field (unless it is rephased by another 
strong pulse). The ability to prepare an ideal inverted 
medium for classical information processing has been in- 
vestigated experimentally (30| . 



VI. PHASE MATCHING 

The treatment so far has involved only one spatial di- 
mension. One way to consider the effect of phase match- 
ing is by extending to a 3D treatment in the paraxial ap- 
proximation. In this case we replace a(z,t) ~f a{z,kt,t) 
and cr_(z,A,t) cr_ (z, p. A, t) where kj = {kx,ky) is 
the transverse wavevector and p = (x, y) is the trans- 
verse position. Our linearized Maxwell Bloch equations 
for the ground state become 



For the excited state Maxwell Bloch, the same procedure 
gives 

^L)^(z,kt,A,i) = iAbl{z,kt,A,t)-id{z,~ktim 
^d{z,kut) = ^ / Dl[z,^kuA,t)dA. (24) 

oz 2.TT 

In the situation where the 7r-pulse is applied off axis the 
phase of the 7r-pulse depends on the transverse position 
leading to the transformation 



D^{z, p, A, t) ^ Dg{z, p, A, t) exp(2^k^.p), (25) 
or after fourier transforming 



D,{z, kt. A, t) <~ Dg{z, kt - 2k^, A, t). 



(26) 



With the RASE pulse sequence described in Fig. [21 the 
phase of the first tt pulse doesn't matter. The atoms are 
all in the ground state before the pulse and assuming a 
perfect tt pulse will end up in the excited state afterward 
regardless. Any small coherent excitation caused by im- 
perfect TT pulses can be ignored, it will quickly dephase 
because of the inhomogeneous broadening and will not 
be rephased as an echo until after the second tt pulse 
which is outside the region of time of interest. The ASE 
caused by the inversion due to this 7r-pulse will be spa- 
tially multimode, with the amount of ASE in a particular 
mode determined by the gain experienced traversing the 
sample. 



Write Process Read Process 



— 13> — 13> 




— 13) — 13> — 13) 




FIG. 4: (Color online) a) DLCZ protocol showing write and 
read process, b) Modified protocol. The inhomogeneous 
broadening of the |1)-|2) transition now leads to an increase 
in bandwidth. 

Suppose we set a detection system to look at the ASE 
produced with wavevector Uase, from Eqs. [23] and 1241 
We can see that the light with this wavevector is entan- 
gled with the atomic excitation with mode — Icase- The 
TT-pulse transfers this to the wavevector —Icase + 2k7r 
according to Eq. [26l Equations [2T] and [22] connect 
atomic and optical modes with the same wavevector so 
we have that the wavevector for the RASE is knASE = 
-kASE + 2k^ or 

kASE + kRASE = 2k^. (27) 

This is the same phase matching condition as a two 
pulse photon echo, kjnput+kocho = 2k^. While this phase 
matching condition is valid outside the paraxial regime, 
the only way to achieve phasematching is with the beams 
colinear or close to conlinear, because the ASE the RASE 
and the tt pulse must all be at the same frequency. 
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to a set of excited state atoms (|1)) spontaneously emit- 
ting into the level |2). The emitted optical field is then 
steered elsewhere for ent ang lement generation with an- 
other ensemble of atoms [13 . Once entanglement is gen- 
erated between two ensembles, one wishes to read out 
one ensembles atomic field to a photonic field in order 
to implement entanglement swapping |19i |. For the read 
process, state |2) becomes the excited state and state |1) 
the ground state. 

One problem with this is the inhomogeneous broaden- 
ing of the |1)-|2) transition causes dephasing limiting the 
time separation between the writing and reading process. 
A modified DLCZ protocol, in close analogy with RASE, 
would overcome this problem. A rephasing pulse on the 
|1)-|2) transition utilizes the inhomogeneous broadening, 
now increasing the bandwidth of the process rather than 
reducing the time separations. The sequence of events for 
this modified DLCZ protocol are shown in Fig. [Hd. It is 
worth noting that the modified DLCZ protocol does not 
have the same issue with echo efficiency as the two level 
scheme because the classical coupling field can be altered 
meaning that the ensemble can be optically thin for the 
writing process and thicker for the reading process. 

The phase matching conditions for the modified DLCZ 
protocol will be the same as given in Eq. [27] for RASE. 
However with a Raman transition it is the wavevector dif- 
ference for the two optical fields that is important. This 
means that one has alot more freedom in the implementa- 
tion because one isn't restricted by the requirement that 
Lu — ck, as one is in the two level case. 



VII. CONCLUSION 



A. DLCZ 

It is interesting to consider the relationship of the cur- 
rent scheme with the DLCZ protocol [!§]. The DLCZ 
protocol involves the creation of entanglement between 
distant ensembles. The relevant energy level diagrams 
are shown in Fig. [4^. Once the level |3) has been adiabat- 
ically eliminated the write process is formally equivalent 



In conclusion we have shown that rephased amplified 
spontaneous emission has strongly non-classical correla- 
tions with the original amplified spontaneous emission in 
the optically thin regime. This leads to the possibility of 
a modified DLCZ protocol, where the problem of dephas- 
ing due to inhomogeneous broadening of the hyperfine 
transitions is solved by a rephasing pulse, increasing the 
bandwidth of the process. 
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